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Abstract 

Every orthonomic system of partial differential equations is known 
to possess a finite number of integrability conditions sufficient to ensure 
the validity of all. Herewith we offer an efficient algorithm to construct 
a sufficient set of integrability conditions free of redundancies. 

1 Introduction 

1. Somewhat surprisingly, existing literature on orthonomic systems lacks 
an effective construction of a provably irredundant sufficient set of integra- 
bility conditions. The aim of this paper is to fill this gap. 

A subproblem attracted much attention in polynomial elimination the- 
ory. When computing Grobner bases through the famous Buchberger al- 
gorithm [4], S'-polynomials arise as analogues of integrability conditions. 
The problem of minimizing the number of ^-polynomials was already ad- 
dressed by Buchberger [5], with many later developments ([35, 8, 6, 1]). 
The basic idea, exploiting syzygies, later migrated back to differential al- 
gebra (Boulier [3]) and Riquier theory (Reid's school [25, 27]). Within the 
syzygy approach one detects superfluous "critical pairs" (^-polynomials or 
integrability conditions) and removes them sequentially. In practical im- 
plementations such as Wittkopf 's dissertation [36] detection of superfluous 
integrability conditions turns out to be nearly exhaustive (see examples at 
the end of the paper) . The solution presented in this paper is of rather dif- 
ferent nature and amounts to a direct construction of the irredundnat set of 
nontrivial integrability conditions, making explicit all remaining freedom of 
choice. 
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2. Present developments of formal integrability theory are, to a great ex- 
tent, driven by computer algebra applications, especially solution of large 
systems of overdetermined PDE connected with computation of symmetries, 
conservation laws, and other invariants of PDE (see surveys [11, 13]). Since 
input systems consisting of hundreds of equations are not uncommon, effi- 
ciency of the algorithms is an important issue. 

Initially (Riquier [24], Janet [14]) the basic question was which coeffi- 
cients of Taylor expansion of a solution could be chosen arbitrarily (para- 
metric derivatives) and which were then uniquely determined by the system 
(principal derivatives). As is well known, hidden dependences between para- 
metric derivatives lead to integrability (or compatibility) conditions. A sys- 
tem with or without unsatisfied integrability conditions is said to be active 
or passive, respectively. The procedure of augmenting an active system with 
its integrability conditions is called the completion. The augmented system 
is not necessarily passive, since new integrability conditions can emerge. 
However, repeated completion is guaranteed to stop after a finite number of 
steps under fairly general assumptions (Tresse [31]; see [22] for an overview). 

Conventional wisdom says that computing integrability conditions amounts 
to taking cross-derivatives. But the notion of cross-derivative applies only 
to orthonomic systems (ones resolved with respect to "highest" derivatives) . 
Moreover, integrability conditions can depend substantially on the way the 
system is resolved (as opposed to the Cartan and Spencer geometric theory 
of involutivity [20] and the recent theory of Mayer brackets [15, 16]). On 
the other hand, if we accept all the unpleasant consequences, as we do in 
the present paper, we find ourselves placed in an environment tailored for 
easy and efficient implementation of reduction. Reduction is a procedure to 
compute a normal form modulo identifications following from the system. 

3. As usual, this paper deals with the infinitely prolonged system 
which consists of individual equations of the input system S differentiated 
with respect to every combination of independent variables. Under a suit- 
able ranking of derivatives, every equation of S°° comes out resolved with 
respect to a principal derivative, which substantially simplifies the procedure 
of reduction with respect to The main technical difficulty is that reduc- 
tion is not unique unless the system £ is passive, which is not guaranteed 
before the completion algorithm is finished. To circumvent this problem, we 
consider a triangular subsystem £' of S°° with the goal to show equivalence 
of £' and 

Our main result says how to locate nonignorable integrability conditions 
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in the monomial ideal(s) generated by the system. Namely, we determine 
substantial integrability conditions at a principal derivative u k from con- 
nected components of certain subset X^ of the monomial ideal ordered by 
divisibility. The idea bears certain surprising resemblance to that of "subcon- 
nectedness" (see [34, 17, 33] and references therein) having roots in Grobner 
basis theory as well, although a closer look reveals substantial differences. 

It is easily seen that u k ^ must be a cross-derivative in order to possess 
nontrivial integrability conditions. This observation immediately leads to 
a canonical and transparent construction of a sufficient set of integrability 
conditions (Construction 4.4). Next, this set is shown to be irredundant 
in a sense that no integrability condition can be ignored. Considering the 
"staircase diagram" associated with the monomial ideal, the main result of 
this paper provides identification of the vertices where nontrivial integrabil- 
ity conditions reside. This opens the door to asking and answering various 
combinatorial questions, which is however out of the scope of the present 
paper. 

It is often argued that orthonomic essentially means linear from the 
practical perspective, since arbitrary nonlinearities can occur at later stages 
of completion. Let us stress that results of this paper rather loosely depend 
on orthonomicity, since for the integrability conditions to show up it is not 
necessary that the system be explicitly resolved. Computation of derivatives 
of implicit functions being little problem, difficulties lie with implementing 
reduction. In case of polynomially nonlinear systems a help comes from the 
theory of triangular systems (see [12] for an overview). That said, we leave 
this issue to a further study. 

An extended abstract of the previous version of this paper appeared in 
Proceedings of the GIFT 2006 conference [18]. Another abridged exposition 
will be made through the book [32]. 

2 Orthonomic systems 

In this section we recall standard facts and fix our notation. We denote 
by II = {u 1 , . . . , u 1 } a set of dependent variables and by X = {x\, . . . ,x n } 
a set of independent variables. Consider the free commutative monoid X* 
over X. An arbitrary element /i 6 X* is of the form = x^ 1 ■ ■ ■ x 7 ^ 1 , r\ > 
0, • • • , r n > 0, and will be called a Janet monomial [14]. A derivative 
d ri+ '" +Tn u k / dxY ■ ■ ■ dx 7 ™ can be identified with a pair (u k , x^ 1 ■ ■ ■ x r ^ ) <G 
U x X*. It will be convenient to denote derivatives as -u^, G X*. Dependent 
variables u k can be idetified with derivatives u k of order 0. 
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Elements of X U (It x X*) bijectively correspond to local coordinates on 
an appropriate infinite-dimensional jet space J°° [2, 20, 28]. For the purpose 
of understanding the present paper it is sufficient to think of J°° as an 
infinite-dimensional space equipped with coordinates indexed by elements 
of X U (IX x X*). Smooth functions are defined as mappings J°° — > R that 
(locally) depend on only a finite number of coordinates. For each x € X, the 
total derivative 

can be viewed as a vector field on J°° (or differentiation of the R-algebra 
of smooth functions on J°°). Observe that D x acts on a derivative by 
multiplying the Janet monomial li by x. As is well known, total derivatives 
commute. For every v £ X* the corresponding composition of total deriva- 
tives is denoted by D v . 

Definition 2.1. We denote by < the relation of divisibility of monomials 
and by li/v the quotient of monomials fi, v € X* whenever li > v. If so, then 
the derivative u£ is said to be lower than the derivative vJ^. We say that v 
is strictly lower than li, writing v < li, if v < [i and v ^ fi. 

Essential in Riquier's theory is a suitable ordering of derivatives com- 
patible with differentiation. 

Definition 2.2. A ranking is a linear ordering ^ of the set U x X* of 
derivatives such that 

P -< D x p, 

p^p' D x p< D x p 
for all p,p' £ U x X* and every x € X. 

Obviously, we have the implication q < p q ■< p, but not the converse. 
It is an easy consequence of Dickson's lemma (see also Thomas [30]) that 11 x 
X* is a well-ordered set, in particular, every decreasing chain is finite. This 
enables transfinite induction through the set UxX*. A complete classification 
of rankings has been obtained by Rust [25, 26]. 

Now we take into consideration a system £ of finitely many partial dif- 
ferential equations, resolved with respect to diverse derivatives 

< = ** (i) 
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where <££j are smooth functions on J°°. By domE we denote the set of all 
derivatives, appearing on the left-hand side of at least one of the equations 
(1) (even though domE is determined by the form (1) rather than by E 
itself). 

The basic object of interest in formal integrability theory is the asso- 
ciated infinitely prolonged system E°° consisting of all possible differential 
consequences: 

u%, = D v **, uGX*. (2) 

Derivatives u k G appearing on the left-hand side of some equation from E°°, 
i.e., belonging to dom E°°, are said to be •principal. The other derivatives are 
said to be parametric. Thus, a derivative is principal if it either belongs to 
domE or is a derivative of some derivative from domE. 

Now, each equation of the input system E is supposed to be resolved with 
respect to the highest rank derivative it contains, and the right-hand sides 
are supposed to be free of principal derivatives. Summing up, the system (1) 
is assumed to be orthonomic in the following sense: 

Definition 2.3. A system of equations E in the form (1) is said to be 

- triangular, if for every derivative q £ dom E there is exactly one equa- 
tion with q appearing on its left-hand side; 

- normal with respect to a ranking ^, or ^-normal, if every derivative 
q £ domE is ^-maximal in its equation; 

- autoreduced, if no principal derivative occurs on the right-hand side of 
any equation; 

- orthonomic with respect to a ranking ^, or ^-orthonomic, if it is 
triangular, ^-normal, and autoreduced. 

The following easy observation is crucial to Riquier's theory: The in- 
finitely prolonged system E°° is normal when E is normal. However, the 
property of being triangular is usually lost in E°° (which is why integra- 
bility conditions occur). Having autoreduced right-hand sides of E will be 
useful in the sequel, while for E°° no such property is needed. 

We came to the point where division between multiplicative and non- 
multiplicative variables [14] enters the discourse in orthodox exposition of 
the Riquier-Janet theory. Formalized by Gerdt and Blinkov [9, 10] (these 
works triggered an active thread of research in polynomial elimination the- 
ory), so-called involutive divisions became a standard tool to prescribe a 
unique right-hand side to every principal derivative from domE 00 . Selecting 
a unique right-hand side is convenient, but not absolutely necessary (Reid's 
school [21, 23] succeeds without it). 
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3 Reduction subsystem 



Here we take a another path to uniqueness. Namely, we consider an arbitrary 
triangular subsystem £' of E°° with dom £' = dom S°° and call it a reduc- 
tion subsystem of £°°, since it provides us with a unique reduction. Then 
completion (see Introduction) may be regarded as a procedure to establish 
equivalence of £' and £°°. All algorithms to follow actually refer only to a 
finite part of hence its infiniteness does not hamper computability. 

Construction 3.1. For each u k ^ G dom£°°, choose arbitrary £ G X* such 
that u*/£ G dom£ (so that exists) and put = £>£<J>^. The trian- 
gular system of equations 

< = «j€domE°° (3) 

is the reduction subsystem sought. 

This construction gives us a considerable freedom of choice, measured by 
the number of elements in dom£ that are lower than v, 1 ^ (see Definition 2.1). 

Given an expression F depending on a finite number of derivatives, one 
can apply equations of the reduction subsystem (3) as substitutions to obtain 
an "equivalent" expression SF without dependence on principal derivatives. 
In each step, the highest rank principal derivative p = u k the expression F 
actually depends on is substituted by the corresponding expression from 
Construction 3.1. Such steps can be repeated while F depends on principal 
derivatives. There can be only a finite number of these steps since -< has the 
descending chain property. Hence the reduction procedure is algorithmic. 
Effective implementations are available, see, e.g., Wittkopf [36]. 

Reduction S is an IR-algebra homomorphism C 00 / 00 — > C 00 / 00 and 
satisfies S o S = S. 

Applying reduction S on the right-hand sides ^ of system (3.1) we ob- 
tain an autoreduced reduction system (see Definition 2.3), whose the right- 
hand sides ^ k depend only on parametric derivatives, hence we have simply 
Su k = The autoreduced system generates the same reduction as the 
unreduced one. Actually, the previous version of this paper (see the extended 
abstract [18]) depended on use of an autoreduced reduction subsystem. How- 
ever, autoreduction is no longer necessary in practical implementations (see 
Remark 6.1). 

Our main result below (Theorem 4.8) shows that reduction and total 
derivatives D x on J°° satisfy 

S o D x o S = S o D x . 
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The next example demonstrates that we have no such property until we 
know that £ is passive. 

Example 3.2. A simple example of an active system £ with SD X SF / 
SD X F is 

Ux = f(u), u y = g{u). 
Let the reduction subsystem £' contain the equation 

rather than its alternative u xy = D y g (the ranking -< can be arbitrary). For 
F = u y we obtain 



(dg_ \ = dg_ 



SD x Su v = SD x q = S ( -—u r = — — f, 



df \ df 

SD X Uy = SU X y = SDyf = S[ 7^2/ J = ^9 '■ 



Observe that SD x Su y = SD x u y is exactly the integrability condition 

df &g 

-£-9 = ^~f 
Oil ou 

for the system £. 



4 Integrability conditions 

Henceforth we fix a reduction subsystem £' of £°° such that domS' = 
dom £°° as in the preceding section. As above, S denotes the reduction with 
respect to £'. Integrability conditions, investigated in this section, measure 
the nonequivalence of various ways of prolongation. 

Definition 4.1. For every principal derivative u k (i.e., u k £ domE°°) we 
introduce the principal subset X^ as the set of all monomials ^ / 1 such 
that <G dom£°°. Thus, elements of X^ are principal derivatives strictly 
lower than u^ (see Definition 2.1). 

Definition 4.2. If in £ there is an equation of the form u k = <3?^ such that 
the principal subset X^ is nonempty and ( G X*, then the condition 

^ = SD^Su^ (4) 
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is called an integrability condition of the first kind at the point vf^. 
For every pair £,77 € X^, the condition 

SDeSufo = SD n Sul h (5) 

is called an integrability condition of the second kind at the point u k . 

Let us remind the reader that £ is not necessarily a subset of £'. This 
explains why integrability conditions of the first kind are needed. 

According to Definition 4.2, integrability conditions at the point u k ^ are 
satisfied if any two possible ways of obtaining the value Su k ^ lead to one and 
the same result. It is well known that all such integrability conditions follow 
from a finite subset. The main goal of this paper is to reestablish this result 
in an irredundant way. 

For every principal derivative u*, consider the principal subset X k or- 
dered by the divisibility relation < . The following definition is quite standard 
and reflects properties of the graph of the ordered principal subset. Let « 
denote the reflexive, symmetric, and transitive closure of the ordering <. In 
other words, p ~ q if and only if there exists a finite sequence of monomi- 
als z r , . . . , Z2s+i G X^ such that p = z±, q = Z2s+i and Z2j-i < £2j whereas 
Z2j > z 2j+i for every j = 1, . . . , s. Since « is an equivalence relation, we have 
got a partition X^/w of the set X^j into equivalence classes [x]~, x G X^. 

Definition 4.3. A pair of elements p, q of the principal subset X^ is said 
to be connected if p ~ g. Equivalence classes with respect to ~ are called 
connected components of X^. The principal subset X^ is said to be connected, 
if it consists of a single connected component, otherwise it is said to be 
disconnected. 

Construction 4.4. For every G dom £ with nonempty principal subset 
X^j choose one integrability condition of the first kind (4) , 

^ = SD^Su^, 

where £ G X^ is arbitrary. For every u k G dom£°° such that the principal 
subset X^ consists of s connected components [£i]~, . . . , [£ s ]~ with s > 1, 
choose arbitrary representatives £1, . . . , £ s of these components and consider 
integrability conditions of the second kind (5) in the form of a chain of 
equations 

The set of integrability conditions obtained in this way will be called a 
sufficient set of integrability conditions. 
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Clearly, each £ of Construction 4.4 can be chosen so that is mini- 
mal in X^, hence belongs to domS so that we can replace SD^Su 1 ^^ with 

SD^K^ to obtain conventional integrability conditions in the sense of the 
following definition. 

Definition 4.5. Assuming u^,^^. G domS, integrability conditions of 
the form 

are said to be conventional. 

Remark 4.6. Obviously, every integrability condition can become conven- 
tional at the cost of enlarging the system £ by appropriate reduced equations 
from 

Our immediate goal now is to show that every sufficient set resulting 
from Construction 4.4 implies all the other integrability conditions of Def- 
inition 4.2. The following lemma is the key. Let var F denote the finite set 
of all variables (independent variables and derivatives) a smooth function F 
depends on. 

Lemma 4.7. Let x be an independent variable, a G X* a monomial, and F a 
function of independent variables and parametric derivatives. Let SD T SD x p = 
SD XT p for every derivative p G varF and every monomial r < a. Then 
SD a SD x F = SD xa F. 

Proof. We have 

D a (FG)= CD p F-D T G, 

pr=a 

for suitable constants c£ T . Applying SD a to 

r) F r)F 
D * F = E SD X F= -g-SD x q, 

(jGvarF qevar F 



we get 



SD a D x F= Y, c - SD p{^)- SD ^ 



gGvar F pr=a 
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whereas 

SD a SD x F= E CSD p ( ^pj- SD T SD x q. 

gGvar F pr=cr ^ y ' 

These two expressions coincide since SD T SD x q = SD TX q holds by assump- 
tion for all q and r < a. 
□ 

Theorem 4.8. Suppose that the reduction subsystem £' (see Sect. 3) sat- 
isfies some sufficient set of integrability conditions as in Construction 4-4- 
Then 

(i) for all u p G dom S°° and all £ < // we have 

Su* = SD^Sn^; (7) 

(ii) all integrability conditions in the sense of Definition 4-2 hold true 
(meaning that £ is passive); 

(iii) for every monomial £ and every smooth function f on J°° we have 

SDtf = SD^Sf; 

(iv) manifolds an d £s' coincide. 

Proof. To prove (i) we proceed by induction with respect to p = u p . If 
exists and is parametric, then (7) is satisfied trivially since Su k ^^ = u p ^. It 
remains to deal with the case when u^,^ exists and is principal, i.e., the case 
of £ <G X^. To start with, we consider p = u k ^ minimal with respect to the 
ordering -<;. Then (7) holds true in a trivial way, since = in that case. 

To perform the induction step, let us consider an arbitrary derivative 
p = u k l assuming validity of (7) for all q -< p. We shall prove 

SD Sul, = SD p Sul /p (8) 

for all a, p G X^. We start with the case of a, p belonging to one connected 
component, i.e., amp. Obviously, this case can be reduced to p < a by the 
definition of ~. But then it can be further reduced to p <\ a meaning that 
p < a and a j p is a variable, since < is the reflexive transitive closure of <3 
in X^. Therefore, let x be an independent variable x such that p = xa. To 
prove that a = p, we establish equalities 

SD a Sul /a = SD a SD x Su k /xa = SD xa Sul /xa = SD p Su k /p . 
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The first equality follows from Su^, = SD x Su^ X(T , which is (7) for p/a -< p 
and £ = x, therefore holds by induction assumption. 

To prove the second equality we apply Lemma 4.7 to F = Su k i ^ xa . Let 
us verify the assumptions. Consider an arbitrary monomial r < a and q € 
varF. Then q -< uj /x<r , whence D XT q ~< D XT u* /xa = u k r/a ± u* a/(J = 
u^. By induction assumption, (7) holds with replaced with D XT q and 
£ with x, meaning that SD XT q = SD x SD T q. Having verified assumptions 
of Lemma 4.7, we have the second equality. The third equality is obvious. 
Thus, (8) holds for arbitrary p, a belonging to one connected component of 

We are left with the case when a, p belong to different components. But 
if the set J of integrability conditions is sufficient in the sense of Construc- 
tion 4.4, as it is supposed to be, then 3 contains an integrability condition 
SD a Su^ a , = SD p Su^ pl with some other p' « p and a' « a from the same 
components, and then (8) holds for a',p' by assumption and then for a,p 
by transitivity. 

This means that we have one and the same value SD„Su k , = SD n Su k , 

At/o" p HIP 

for all a, p € X^. To establish (7), it remains to show that this common value 
is also equal to Su^. If domE, then Su* = SD^Su 1 ^^ for some ^elj 
by construction of the reduction system E' (Construction 3.1). If € dom E 
and X£ = 0, then (7) is void. If u£ € domE and X* ^ and Su* = 
then the sufficient system involves an integrability condition of the first kind 
<I>^ = SD^SuK^ for some £ E X^. Finally, if <G domE and 7^ and 

Su* 7^ 3>k, then (7) follows from the same Construction 3.1 again. Thus, 
statement (i) is proved. 

Statement (ii) follows immediately from (i) or (8). Statement (iii) holds 
for all functions / if it holds for all derivatives u k , and then it follows from (i). 
Finally, by (ii) every equation from system E°° becomes an identity when 
reduced with respect to S. This means that systems E°° and E' are equiva- 
lent (follow one from another). Hence statement (iv). 
□ 

Remark 4.9. Operators SD X S on the full jet space J°° are K-linear and 
satisfy the Leibniz rule, hence they are vector fields. Moreover, they com- 
mute, since [SD x S,SD y S]f = SD x SD y Sf - SD y SD x Sf = SD x D y Sf - 
SD y D x Sf = by Theorem 4.8(iii). Hence the full jet space J°° equipped 
with vector fields SD X S, x G X, is a diffiety in the sense of [2]. Now, by 
Theorem 4.8(iv) we have C°°E s °o = C°°£ s / ^ C°° J 00 / Ker S ^ SC^J 00 . 
Hence SD X S induce well defined operators on manifold £s«> as well, turning 
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it into a diffiety. 



In an attempt to convey the sense of the method we conclude this sec- 
tion with a simple example in dimension two. For less trivial examples see 
Section 8. 

Example 4.10. Consider the following system E: 

Uxyyyy = C> U X xyyy = f i U X xxyy = 9i U X xxxy = h, 

where e, /, g, h are arbitrary functions of parametric derivatives. Figure 1 
shows the ordered set dom S°° of principal derivatives placed within a coordi- 
nate system. Symbols e, /, g, h denote the four generating derivatives u xyyyy , 
^xxyyyi ^xxxyyj ^xxxxy ^ doinE, respectively. The bold dot at u xxxxyyyy de- 
notes a typical principal derivative. Thick lines show the principal subset 
Xxxxxyyyyi which is obviously connected. Actually, one easily sees that all 
principal subsets are connected except 



T 

•^xxxxyy 


= {u, 


X 

•^xxxyyy 


= W. 


X 

•^xxyyyy 


= W. 



which consist of two isolated points each. Correspondingly, each of the 
derivatives u xxxxyy , u xxxyyy , u xxyyyy harbours one nontrivial integrability 
condition. They are, respectively, 

SD x e = SD y f, SD x f = SD y g, and SD x g = SD y h. 



5 Cross-derivatives 

In this section we find an alternative description of the sufficient set suitable 
for effective implementation. 

Construction 4.4 leaves us with the problem to find all principal deriva- 
tives with disconnected principal subset. Indeed, a nontrivial integrability 
condition of the second kind at a point n^j exists if and only if there are at 
least two distinct connected components in . To extend this line of argu- 
ment further, let us consider different possible descriptions of the quotient 
sets X£/«. 

Let B be a subset of X^j such that every element of X^ is connected to an 
element of B in the sense of Definition 4.3. Then, obviously, every connected 
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Figure 1: A typical principal subset (Example 4.10) 




component intersects with B. In particular, the quotient set X^/« is the 
same as the quotient set B/^b, where is the equivalence relation on B 
inherited from the relation ~ on X k There are two natural choices for B, 
which lead to two alternative descriptions of X^/«: 

(a) minX^ = the subset of minimal elements in X^; 

(b) maxX^j = the subset of maximal elements in X^. 

Let _/V"( fc ) denote the set of all monomials fi such that u k G domS. Assum- 
ing _/V"( fc ) ordered by divisibility, let = minN^ denote the set of all 
minimal elements in N^ k \ Obviously, the subset minX^ coincides with the 
intersection X^ n miniV"( fe ) = X^ n M^ k \ Define a reflexive and symmetric 
relation j on minX^ by p f q if lcm(p, q) G X^, i.e., if lcm (p, q) is a proper 
divisor of u k . . 

Elements of maxXj are quotients \ijx with x G X an independent vari- 
able such that the derivative u^, is principal. To simplify reasoning, we 
identify maxX^j with a subset of X. Define a reflexive and symmetric re- 
lation I on maxX^ C X by x [ y if there exists a G (equivalently, 
a G Af( fc )) and a < fi/x, /i/y. The same relation J. can be defined by x j y if 
x = y or the derivative Up/ xy is principal as well. 

We have the following obvious lemma. 
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Lemma 5.1. The inherited equivalence relation ~ min x fc coincides with the 
transitive closure |* of |. TTie inherited equivalence relation ~ max x fe coin- 
cifies mi/t the transitive closure [* of [. 

Corollary 5.2. We have bijections 

Xj/ra <-► minX^/T* <-> maxX^/j*. 

Proposition 5.3. Xei 6e a principal derivative such that the principal 
subset X^ contains nonequivalent elements a ft r. Then \i = lcm(<r, r). 
Elements <t,t can be chosen lying in minX^. 

Proof. Since a, r G X^j, we have lcm(<r, r) < /U. If lcm(<j, r) < fi, then obvi- 
ously a ~ r, contradicting the assumptions. The last statement follows from 
the fact that every connected component intersects with minX^. □ 

Now we introduce cross-derivatives as the "least common derivatives." 

Definition 5.4. A cross-derivative is a derivative "U\ cm ^ aT y where u^,u^. G 
dom S and a, r do not divide one another. 

By Proposition 5.3, nontrivial integrability conditions of the second kind 
can be found only at cross derivatives. Hence the well-known result that the 
number of such integrability conditions is always finite and less or equal to 
\p(j> — 1), where p is the number of equations in the system S. 

Of course, a cross-derivative gives rise to integrability conditions (6) if 
and only if it satisfies the following nontriviality condition: 

Definition 5.5. A cross-derivative is said to be trivial if the principal 
subset X?, is connected. Otherwise it is called nontrivial. 

Example 5.6. Generalizing Example 4.10, consider an arbitrary system of 
r equations in two dimensions such that domS consists of incomparable 
derivatives (with respect to < of definition 2.1). It is an easy exercise to 
show that of the r(r — l)/2 cross-derivatives only r — 1 are nontrivial. 

Let us finish this section with some remarks concerning visualization of 
the relation j. The monoid X* can be visualized as the n-dimensional grid 
N n C R n , where N = {0,1,2,...}, via the correspondence x r ^ ■ ■ ■ x r n n «-> 
(n, . . . , r n ). Given a monomial /i € X* the cone generated by \i is defined to 
be C(fx) = {/j,u \ v G X*}. A union of cones in X* is called a monomial ideal 
(see, e.g., [19]). For each k, we have 

{a <E X* | u k a G dom£°°} = J C(ji). 

Gdom E 
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Hence, to every infinitely prolonged system S°° there corresponds a collec- 
tion of monomial ideals, one for each k, consisting of principal derivatives 
u k . with one and the same k. 

Monomial ideals are usually visualized by staircase diagrams in R n . In 
W l , every point (zi,..., z n ) € N™ generates the corner 

C(zi, ...,z n ) = {{xi, . . . ,x n ) e R n | Zi < Xi for all i}. 

A union of corners, which is an unbounded orthogonal (usually non-convex) 
polytope with vertices in integer points N n C M n , is called a staircase dia- 
gram. An oriented edge between two integer points p = (z±, . . . , Zj, . . . , z n ) 
and q = (z±, . . . , zi + 1, . . . , z n ) is called a direction from p to q. A square 
bounded by four adjacent edges is called a tile. An xy-tile is a tile parallel 
to the xy-plane. Two of the bounding directions end in a common point, 
which will be called the vertex of the tile. Now, on the staircase diagram 
lies in, maxlj can be seen as the set of all directions that lead to p,. Two 
distinct directions x, y G maxX then satisfy x [ y if and only if the staircase 
diagram contains the xy-tile with the vertex /i. 

6 The algorithm 

Before proceeding to more substantial examples, let us finally present the 
procedure to find a sufficient set of integrability conditions. Below the sym- 
bol # denotes the number of elements in a finite set and var(p) = {x G X | 
x divides p} for p <G X* is the set of all variables to occur in a monomial. For 
clarity, we present two separate algorithms, one for integrability conditions 
of each kind. Both algorithms share partitioning of the system X into sub- 
systems such that dom contains derivatives u 1 ^ of u k . Sets domE^'' 
being denoted by N^ k \ their minimal elements are then collected in subsets 
M( fe ) C AT«. Of course, the algorithms can share the k loop. 

Algorithm 1 to compute integrability conditions of the first kind is very 
simple. To each non-minimal element u k £ \ there corresponds 
exactly one integrability condition of the first kind. 

Before explaining Algorithm 2 to compute integrability conditions of the 
second kind, let us consider various descriptions offered by Corollary 5.2. 
There is no upper bound for the size # minX^ since the number of equations 
in the system £ can be arbitrary, whereas ^ max X^j is bounded by the 
number of independent variables #X. During completion, # min X^ grows as 
new equations are added to the system, while # max 1 k remains essentially 
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Algorithm 1 Integr ability conditions of the first kind 
Input: S. 

Output: the set IC^ of integrability conditions of the first kind. 

1 IC4 1} := 

2 for all k do 

3 AT(fc) := {„ G X* | u k v € dom £} 

4 M^ fc ) := the set of minimal elements in with respect to < 

5 for jU € JV ( fc ) \ mW do 

6 select arbitrary v £ such that v < n 

7 adjoin expression — D^/^* to IC^ 

8 end for 

9 end for 

10 return IC^ 



stable. This is why Algorithm 2 uses the relation J. on maxlj rather than 
the relation f on minX^. 

Now, Algorithm 2 works as follows. On line 4 cross-derivatives are com- 
puted, using only minimal derivatives collected in M^ k \ This is justified by 

(2) 

Proposition 5.3 as trivial cross-derivatives do not contribute to IC^ . 

To explain lines 6 and 7, observe that the reflexive and symmetric re- 
lation { is conveniently represented as the union of relations [ a , where a 
runs through and | CT is defined by x [ a y if x = y or a < fi/x,[i/y. 

Obviously, each is an equivalence relation. Moreover, the associated par- 
tition of X C maxXj has only one nontrivial class, apart from one-element 
sets, and this class can be identified with var(///<r). Hence, the partition 91 
corresponding to the transitive hull [* is the least partition such that every 
var(///cr) is a subset of some class of 91. Computing such a least partition 
amounts to joining all incident subsets. Alternatively, 91 can be described as 
the set of connected components of the hypergraph (maxX^, 971), where the 
set of hyperedges is 971 = {var(^/a) | a in M^} and the set of vertices is 
max X^ = |J 971. Known algorithms are capable of labeling connected com- 
ponents of a (hyper)graph in expected time linear in the number of vertices, 
see [29]. The time complexity of Algorithm 2 is estimated in Remark 8.3 
below. An obvious modification is to interlace lines 7 and 8 and break as 
soon as 91 has only one connected component. This will further reduce the 
average running time. 

Of course, time spent on redundancy elimination is secondary in com- 
parison with large variance resulting from different completion strategies. 
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Algorithm 2 Integrability conditions of the second kind 

(2) 

Input: £. Output: the set IC^ of integrability conditions of the second kind 

1 icg } := 

2 for all k do 

3 let , AfW be those of Algorithm 1 

4 := {lcm(z/i, 1^2) | {^1,^2} Q M( fc ) a two-element subset} 

5 for ,u € do 

6 Wl = {var(p/a) \ a G and a < fi} 

7 91 := the set of connected components of the hypergraph (|J Wl, 9JT) 

8 if #01 > 1 then 

9 select arbitrary a±, . . . ,a s G such that var(^/crj) is a subset 
of the ith connected component, 

10 adjoin - D„ /<T , - D m/(T1 to icg } 

11 end if 

12 end for 

13 end for 

14 return ICL 2) 



Vital for selecting the proper completion strategy is knowing the available 
freedom of choice. It is, however, clear that the arbitrary selections made 
on line 6 of Algorithm 1 and lines 10, 11 of Algorithm 2 exhaust the entire 
freedom of choice of conventional integrability conditions (Definition 4.5) 
relative to Remark 4.6. 

Let us discuss the whole completion algorithm now. If all the integrabil- 
ity conditions of Theorem 4.8 are satisfied (i.e., if they reduce to identities 
as explained in Sect. 3), then the system is passive and no further steps are 
needed. Otherwise let £ denote the extended system obtained by resolving 
the non-identical integrability conditions with respect to maximal deriva- 
tives (under the same ranking) and adjoining them to S. The system E is 
afterwords subject to the same procedure of selecting a reduction subsystem 
E' of the infinite prolongation E°°, identifying the integrability conditions, 
etc. Obviously, we then have £ S ' != ^s' an d SoS = S = SoS. 

Of course, the new integrability conditions resulting from extending E 
to E can be of both first and second kind. At the same time yet unresolved 
integrability conditions of the second kind can trivialize, which is very easy 
to detect with Algorithm 2. This explains why proceeding incrementally 
(resolving one integrability condition at a time) lowers the overall cost of 
completion. Moreover, the ratio of trivialized and added integrability con- 
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ditions can be computed quickly and reliably if knowing the new leading 
derivative alone. 

Of course, being a part of a completion algorithm, the implementation 
should keep track of the integrability conditions already satisfied in previ- 
ous steps. But do the once satisfied integrability conditions of E continue 
to hold under the new reduction S? Recall the inductive proof of (7) (the 
core statement of Theorem 4.8 the others follow from) for the system E . As- 
sume (7s), i.e., Su k ^ = SD^Su 1 ^^, at while for all u l v -< u k assume (7g), 

i.e., Su l v = S D^Su l u ^. By (7g) we have SD^SF = SD^F for arbitrary func- 
tion F of derivatives that precede u^, in particular, for F = SuK^. But then 
SD^Su^ = SD^SSu^ = 5£> € SuJ /€ = SSD^Su^ = SSu* = Sv*. 

Thus, we have (7g) at u k as well. 

An important question is whether the completion algorithm eventually 
stops. An affirmative answer easily follows from the Dickson lemma, since 
new integrability conditions can reside only at points u k outside the mono- 
mial ideals generated by domS^. 

We finish this section with a remark on autoreduction. 

Remark 6.1. Certain grounds exist for maintaining the reduction sub- 
system non-autor educed. For example, let the input system E contain two 
equations = Fu, u v = Gu, where F, G are linear differential operators 
with constant coefficients. Let, moreover, /x, v be relatively prime. Then the 
corresponding integrability condition of the second kind at u^ u is nontriv- 
ial, yet automatically satisfied: D u Fu = FGu = GFu = D^Gu on 6^. 
(In polynomial elimination theory this case is covered by the so-called first 
Buchberger criterion [5].) Now, the crucial identity FGu = GFu (which 
only appears in expanded form) is much easier to check before applying any 
reductions for u a <G var FGu = var GFu. To a lesser extent this is so even if 
F, G have non-constant coefficients etc. 

7 Irredundancy 

In this section we prove that Construction 4.4 produces no redundant inte- 
grability condition. By a redundant condition one usually means one that 
can be safely omitted from the checklist, since it is satisfied automatically 
whenever all the others are. To put it more formally, observe that essentials 
of Construction 4.4 depend only on the set P = dom E of derivatives on 
the left-hand side of the input system (1), while functions 3>| on the right- 
hand side play the role of parameters. Let us therefore consider the whole 
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class Sp of orthonomic systems £ with fixed set P = domE, parametrized 
by arbitrary functions $| subject only to requirements of orthonomicity. Ob- 
viously, Construction 4.4 provides a set of integrability conditions applicable 
to every member of the class Sp. 

Definition 7.1. Consider the class (Sp of orthonomic systems £ with a 
fixed set P = domE. Let J be a set of integrability conditions of (Sp. An 
integrability condition I G J is said to be redundant if it is satisfied for every 
choice of right-hand sides $| for which all the other integrability conditions 
J \ {/} are satisfied. The set J is said to be irredundant if it contains no 
redundant integrability condition. 

A chain (6) is to be considered as a sequence of s — 1 integrability con- 
ditions, so that each equality sign determines a separate integrability con- 
dition. 

Remark 7.2. Definition 7.1 implicitly refers to some functional space S to 
choose the right-hand sides $| from. Proof of Proposition 7.3 below only 
requires that S contains all polynomials in the independent variables. 

Proposition 7.3. The set 3 of integrability conditions resulting from Con- 
struction 4-4 is irredundant. 

Proof. To start with, we assume that all integrability conditions from J 
are conventional (see Definition 4.5). Let = SD^/^SQ* G J be such an 
integrability condition of the first kind. By assigning Q k = 1 and & a = for 
all u l a € domS \ we obtain an orthonomic system obviously satisfying 
all integrability conditions except $^ = 1^0 = SD^^SQ^. 

Similarly, consider an arbitrary conventional integrability condition of 
the second kind from J, say 

SD./^l = SD^\ 2 =■■■= SD^l (9) 

at u k ^ E domS°°. Let [£i], . . . , [£ s ] denote the corresponding equivalence 
classes in minX^. Let 1 < r < s be an arbitrary integer and / denote the 
rth integrability condition in the chain, i.e., SD^/^Q^ = SD^/^^Q* . 
By Construction 4.4, 3 contains no more than one integrability condition of 
the first kind of the form = SD^j^^. If such a a exists, let S denote 
[£i] U • • • U [£ r ] or [Cr+i] U • • • U [£ s ] whichever contains a. If no such a exists, 
then let S be one (arbitrarily chosen) of these two sets. 
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For every monomial a = x" 1 • • • x" n G X* we introduce the function of 
independent variables 



X ®1 . . . rjf-l 

■Fai%l> • • • j En) = 



a±\ ■ ■ -a n \ ' 



which obviously satisfies 

DrFa = lF*/T ifr<a, (1Q) 
1 otherwise. 

Turning back to our proof, for every u l a G dom S we assign <& l a according to 
the following simple rule: 

$ l =I F t*/<r if I = k and a G {/x} U S, /j-q 
lO otherwise. 

Consider an arbitrary integrability condition = SD u / a ^, a < u, of 
the first kind from J. The only possibility how the left-hand side can be 
nonzero is when 

(A) a G {fi} U S, !/</*, 

and then it equals D v j a Q* = F) v j a F ii j a = F^/ u . The only possibility how 
the right-hand side can be nonzero is when 

(B) v G {/i} U 3, 

and then it equals the same F^/y. It remains to be checked that conditions 
(A) and (B) are equivalent. Before that we observe that the inequality a < v 
implies 

(C) if a, v < fi, then a, v both or neither lie in S. 

Indeed, under these conditions we have a « v in X^. 

Let (A) be true. The case v = fj, being trivial, consider v < fx. Then also 
a < [i and therefore uEH. But then v G S by (C), giving (B). Conversely, 
let (B) be true. If z/ = then <j G E since S was chosen that way. Otherwise 
v G S, but then cr G E by (C) again. Therefore, (A) is true. Thus, we have 
proved the equivalence (A) 44> (B) and hence validity of all integrability 
conditions of the first kind. 

Now consider an integrability condition of the second kind from J, at 
some u l u . In case of u l v = the integrability condition is (9). However, 
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SD^/£&£ equals 1 if £ £ £ and otherwise, hence all equalities (9) hold ex- 
cept /. Thus we are left with the case of an integrability condition SD v j ax = 
••• = SD u / at <&£ t from J, at some u l u / u k . By (11) and (10), the values 
SD v j cr .^ k ai to be compared at u l u are all zero except when 

(Aj) I = k, ai<u < fj,, Oi G S, 

and then they are D u / a .^. = D^^.F^^. = F^j v independently of i. Let us 
show that conditions (Aj) are mutually equivalent. However, if one of (Aj) 
holds, then / = k and v < //, hence v < fi (otherwise = u^). We have 
<7j < ^ for all j by Construction 4.4, hence o\ ~ ■ ■ ■ ~ at in (although 
not in X($). Therefore all try belong to S. Equivalence of conditions (Aj) is 
thereby established. 

Thus we have proved the proposition in case of conventional integrability 
conditions. But since all assigned during the proof were functions of in- 
dependent variables only, we have simply SD u <& l a = D u <& l a for the reduction 
S of any principal derivative. This means that every integrability condi- 
tion can be identified with a conventional integrability condition. Hence the 
proposition holds for general integrability conditions as well. □ 

It remains to compare our definition of redundancy with that used by 
other authors, notably Rust [25, 27]. Consider the free abelian algebra A k 
over the set of abstract generators of the form D^^g. The total derivatives 
D x act upon the generators, hence upon the whole algebra, in a natural way. 
An integrability condition can be viewed as a difference D^$| — D^^Q* <G 
A k of two generators. Given a finite set J C A k of integrability conditions, 
another integrability condition / = — D^G** is said to be syzygy 

redundant if monomials m < //, integrability conditions Ii = D^./^.^. — 
^v-i/m^m ^ an d integers Cj6Z exist such that 

1 = E ^^w* = E «( d m*1 - (12) 

holds in A k . 

It is clear that if / is syzygy redundant, then it is also redundant in the 
sense of Definition 7.1 for all choices of S (see Remark 7.2). Hence the suffi- 
cient set of integrability conditions resulting from Construction 4.4, proved 
to be irredundant when S contains all polynomials in independent variables, 
is also syzygy irredundant. 
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8 Examples 



The first two examples compare our algorithms to Algorithm 9 from Wit- 
tkopf's dissertation [36]. Wittkopf's algorithm removes syzygy redundancy 
nearly completely. Experiments with randomly generated monomial ideals 
showed that Wittkopf's algorithm can miss r redundant integrability condi- 
tions in case of ideals with 4r generators, but such instances are rather rare. 
A surprise was that Wittkopf's algorithm could be substantially slower. 

Example 8.1. Consider a system X of the form 

Uxyz = /l; U X xz = hi ^"yyz = /3; U X xyy = fi- 

We summarize the work of Algorithm 2 in a table: 







var(^/cr), a = 


max X M /« 


I C (2) 




xyz 


2 2 2 2 

x z y z x y 






x 2 yz 


x 


y 


{x},{y} 


Dxfi = D y f 2 


xy 2 z 


y 


X 


{x},{y} 


Dyh = Dxh 


2 2 

x y z 


xy 


y x z 


{x,y},U} 


Dxyfl = D z f4 



The first column lists all possible cross-derivatives \i. Columns 2-5 corre- 
spond to the four derivatives u c from domS. These four columns list vari- 
ables the monomial fi/a depends on whenever a divides fi, and contain an 
empty space when [ija is not a monomial. The sixth column contains the 
least partition of maxX M generated by the sets occurring in columns 2-5 
(maxX M is the union of these sets). By results of Section 5 this partition 
corresponds to the equivalence relation « on maxX M inherited from X M . Al- 
gorithm 2 also says how to choose the integrability conditions (we omit the 
reduction symbol). In the first and second row the only possibility is that 
given in the last column (when \i = x 2 yz or [i = xy 2 z, each connected com- 
ponent of X M contains a single a). Contrary to that, one of the connected 
components of X XX yyz contains three monomials a, namely xyz, x 2 z, y 2 z. 
Hence, apart from D xy fi = D z fi shown in the table, there are two other 
equivalent ways to write the third integrability condition: D yy h = D z f4 
and D xx f 3 = D z f A . 

The example can be visualized (see end of Sect. 5) as follows: 
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The reader may wish to locate the tile that induces the equivalence relation 
x Tti y at n = x 2 y 2 z. 

Example 8.1 is one of the simplest where Algorithm 9 from Wittkopf's 
dissertation [36] misses one redundant integrability condition. However, Wit- 
tkopf's algorithm depends on a choice of what is called compatible rank- 
ing [36, Def. 10] of syzygies, which itself depends on a choice of a ranking 
for X (which we fix to be x -< y -< z) and a permutation of the set S. In 
Example 8.1, Wittkopf's algorithm has a very favorable ratio y| °f correct 
answers in the set of all 4! = 24 permutations of S. This ratio can be less 
favorable in other examples. 

Example 8.2. Consider a system £ of the form 

^xxy fli "U>xxz jli ^xyy fzi 

Uxzz = fii ^yyz = fbi Uy ZZ = /g. 
We summarize the work of Algorithm 2 in a table: 





x 2 y 


var(p/a), a 

2 2 2 

x z xy y z 


xz 2 


yz 2 


maxX M /« 

I C (2) 






x 2 yz 


z 


V 








{y},{z} 


D z h 


= Dyh 


xy 2 z 




z 


X 






{x},{z} 


Dzh 


= D x f 4 


xyz 2 








y 


x 


{x},{y} 


Dyh 


= D x f§ 


2 2 

x y 


y 


X 








{x},{y} 


Dyh 


= D x f 3 


2 2 

X Z 




z 




x 




{x},{z} 




= D x f 5 


2 2 

y z 






z 




y 


{y},{z} 


Dzh 


= Dyh 


2 2 

xyz 


yz 


y xz 


X 






{x,y,z} 






2 2 

x yz 


z 


yz 




xy 


X 


{x,y,z} 






2 2 




z 


xz 


y 


xy 


{x,y,z} 
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For explanation of the table see Example 8.1. The last column shows the 
unique integrability condition in each of the first six rows and none in the 
remaining three. The corresponding diagram is 




Wittkopf 's Algorithm 9 gives an incorrect number of integrability condi- 
tions (seven) in 540 of the full number 6! = 720 of permutations of S. Thus, 
the ratio of correct answers is only 1/4 now. 

Remark 8.3. In randomly generated examples, Wittkopf 's Algorithm 9 ran 
substantially longer than ours on the same data. Both algorithms take ad- 
vantage of the partitioning S = |J Y,( k \ Let us therefore attempt comparison 
in case of one dependent variable (so that there is no k loop). 

The outer loop of Wittkopf 's algorithm runs over the syzygy system S, 
which has 0(r 2 ) elements, where r = At each run, subset §' C § of 
already executed (accepted or rejected) syzygies is incremented. Processing 
elements s' € S' in Step 3.1 costs time units. Processing pairs of elements 
s', s" G §' in Step 3.2 costs at least #§' time units again (s' and s" are not 
independent). This suggests running time at least 0(r 4 ). 

In our Algorithm 2, as well as have 0(r) elements. The main 
loop 5-12 runs over L^ k \ which has 0(r 2 ) elements. At each run, build- 
ing 5DT on line 6 requires time proportional to ftM^ . Obtaining connected 
components of the hypergraph 9JT on line 7 requires time proportional to 
#MW + #X, where typically #MW > #X. This suggests 0(r 3 ) running 
time. 

Example 8.4. Since the ring M.[d/dxi, . . . , d/dx n ] of linear differential oper- 
ators with real coefficients is isomorphic to the polynomial ring M.[xi, . . . , x n ], 
PDE algorithms can be applied to polynomial ideals as well, in which case 
they compute Grobner bases. 

Caboara et al. [6] used the syzygy approach to minimize the number 
of S'-polynomials when computing a Grobner basis of a polynomial ideal. 
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When their Example 26 is translated into the PDE language, the following 
Janet monomials correspond to maximal derivatives: 

°1 = x 2 x 3 x 4 x ti °3 = x l x 2 x 3i 

O2 = xfx2 x A x i, a i = x i x t x 5- 

As summarized in the following table, Algorithm 2 reveals rather immedi- 
ately that one of the four existing cross-derivatives is trivial: 

fi subsets C M (cr), u = X M /« 

o\ a 2 &3 O4 

x\x\x%X^x\ X\ X2X3X5 X4X5 X2X4X5 {xi}, {X2, x 3, x 4, 25} 

x l x 2 x 3 x 4 x i X 2 X 3 X 4 X 5 X 2 X 4 { x 2, x 3, x 4, x $} 

xfx2X^X4xl 23 X2X4 {X2, 24}, {xs} 

x\x 2 x\x\ x 5 x 2 {x 2 }, {x 5 } 

Thus, we arrive at exactly three nontrivial integrability conditions in full 
accordance with the result of [6, Ex. 26]. 

Finally, we give an example where integrability conditions ordered by 
divisibility form a chain. 

Example 8.5. Let n > 2 be arbitrary. Consider the following n Janet 
monomials in n variables: 



22X3X4- • • x n , 

X^X3X4 - • • X n , 
2 2 

X \ X 2 X ^' ' ' X Tl} 



2 2 2 2 
■ L 1 X 2 X 3 x n-l- 



The cross-derivatives 



XiX 2 X 3 X 4 - --Xn, 
2 2 

x ~y x 2 x 3 x 4' ' ' x ni 

2 2 2 

12 3 4 "^n ? 



222 2 

X 1 X 2 X 3 x n— l^m 



are all nontrivial and form a chain of length n — 1 . 
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9 The poset of nontrivial cross-derivatives 



In this section we explore simplest relations between divisibility properties of 
cross-derivatives and triviality of the corresponding integrability conditions 
of the second kind. Results concerning different variables u k being totally 
independent, it suffices to consider the case of a single variable u and omit 
the upper index k. 

Given an orthonomic system S, let P^ denote the poset of nontrivial 
cross-derivatives in the sense of Definition 5.5 under ordering by divisibility. 
The following proposition shows that minimal elements of Ps coincide with 
minimal cross derivatives of S. 

Proposition 9.1. Let \i be a minimal cross-derivative. Then is nontrivial, 
meaning that jjl G Ps- 

Proof. We remove nonminimal elements of dom X first to ensure that dom £ 
is an antichain without affecting the minimal cross-derivatives. Let \i = 
lcm(a, (3) with a ^ (3 in domS. Supposing that n is trivial, we arrive at 
a contradiction. By triviality of [i we have a ~ (3 in X^. Hence there exist 
&)■■■> C2s+i G X^ such that a < Ci, Ci > C2, C2 < C3, • • • , C2S-1 > C2s, 
C2s < C2S+1, C2S+1 > /3- Consequently there exist a\,...,a s G domS such 
that (21-1, (21, (21+1 G C(o-j). Then lcm(a, a\) < Ci < M, which contradicts 
minimality of /U unless a, o~\ are comparable. But dom S is an antichain, 
hence a = <7i. By repeating the same argument we get a = u\ = 02 =■ ■ • = 
a s = (3, contradicting the assumptions on a, (3. □ 

Hence minimal elements of Ps coincide with the minimal integrability 
conditions in Reid's sense [21]. It follows that nonminimal and nontrivial 
integrability conditions in our sense can be interpreted as a solution to Reid's 
problem of finding an irredundant set of supplementary conditions. 

Turning back to the poset Ps, we show below that it can be an antichain 
(the generic case; all points are minimal) as well as a chain of length less 
than the number n of independent variables (a single minimal point). 

Lemma 9.2. Let v < /i be two nontrivial cross-derivatives. Then X u C X M 
lies entirely within an equivalence class of « at u^. 

Proof. Obviously, X u C X^. Let £,7/ G X u be arbitrary. Since G X M 

and £ < is, rj < is, we have £ w r/ at and the lemma follows. □ 

It easily follows that Ps does not contain a chain of length n, where n is 
the number of independent variables. Yet Example 8.5 shows that Ps can 
be a single chain of length n — 1. 
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In case of general position, P^ is an antichain and all nontrivial cross- 
derivatives are minimal: 

Proposition 9.3. If the generating points of domS C N n have no coordi- 
nate in common, then P^ is an antichain. 

Proof Let a = x^ 1 ■ ■ ■ x^" , (3 = x^ 1 ■ ■ ■ x h £ be such that cij / 6« for all i = 
1, . . . ,n. Let fj, = lcm (a, (3) be nontrivial. However, var(/i/a) and var(/Lt//3) 
are two disjoint sets whose union is the whole X, hence they constitute 
the partition X/~. By Lemma 9.2, there is no room for a nontrivial cross- 
derivative above (i. 
□ 

Interesting combinatorial questions about Ps had to be left aside. Know- 
ing the number of elements of Ps can be useful, the more so since there seem 
to be no similar results under the syzygy approach. For #Pe we have the 
following exact upper bounds in low dimensions: — 1 when n = 2 (obvi- 
ous) and 3 #S — 6 when n = 3 (the same as the number of edges in a planar 
graph). However, computer experiments show that "average" numbers are 
substantially lower. Combinatorial aspects of Pe are also a subject of the 
recent work [7]. 
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